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Aims. To systematically stress a rotationally symmetric 3D magnetic null point by adverting the opposite footpoints of the spine axis 
in opposite directions. This stress eventually concentrates in the vicinity of the null point forming a local current sheet through which 
magnetic reconnection takes place. The aim is to look for a steady state evolution of the current sheet dynamics which may provide 
scaling relations for various characteristic parameters of the system. 

Methods. The evolution is followed by solving numerically the non-ideal MHD equations in a Cartesian domain. The null point is 
embedded in an initially constant density and temperature plasma. 

Results. It is shown that a quasi-steady reconnection process can be set up at a 3D null by continuous shear driving. It appears that 
a true steady state in unlikely to be realised as the current layer tends to grow until restricted by the geometry of the computational 
domain and imposed driving profile. However, ratios between characteristic quantities clearly settle after some time to stable values 
- so that the evolution is quasi-steady. The experiments show a number of scaling relations, but they do not provide a clear consensus 
for extending to lower magnetic resistivity or faster driving velocities. More investigations are needed to fully clarify the properties 
of current sheets at magnetic null points. 
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1. Introduction 

The solar photosphere is threaded by a complicated distribu- 
tion of magnetic field concentrations of opposite polarities. 
Above the photosphere, in the equatorial region, the positive 
and negative flux concentrations typically connect to one an- 
other in a non-trivial way, creating complicated topological pat- 
terns. The simplest approach to investigate the topology of the 
magnetic field is to use magnetogram data as boundary condi- 
tions for potential or linear force free magnetic field extrapola- 
tions. This has been done by a number of authors (IClose et alJ 
I20031 ISchriiver"^ Titld I2003t lLongcope et ail [2003) and is 
found to lead to large numbers of separatrix surfaces divid- 
ing three-dimensional (3D) space into regions of different mag- 
netic connectivity. The analysis show that 3D nulls are present 
in abundance at low altitudes, becoming increasingly rare high 
above the photosphere, and the number distribution depends 
on the distribution patterns of the photospheric magnetic field 
( lLongcope & P arnell 2009). From observations there are also a 
few cases in which 3D nulls have been inferred from the struc- 
ture of the emission in the corona (Filippov 1999), but using this 
indirect evidence for a 3D null is problematic without making 
attempts to derive the underlying magnetic field structure. This 
was recently done by [Masson et al. (2009) who verified their 
magnetic field extrapolations by comparing with the observed 
emission patterns. These independent approaches all show that 
3D nulls are present in the solar atmosphere, and we must there- 
fore investigate their role in the energy release process that takes 
place there. 

Magnetic reconnection is though to play a leading role in a 
large number of different types of magnetic energy release pro- 
cesses where fast conversion of magnetic energy into bulk mo- 
tion, plasma heating and particle acceleration is observed. These 



processes have been observed in a wide variety of situations, for 
example in "explosive" events taking place in the outer solar at- 
mosphere. To understand these phenomena in their various phys- 
ical environments, a detailed understanding of the reconnection 
process is required. For many years two-dimensional (2D) re- 
connection models have formed the basis for this understanding. 
Analytical investigations have shown that various configurations 
with different characteristics are possible, and the driving profile 
controls these characteristic properties to a large extent. One key 
property that can be estimated is the scaling behaviour of the 2D 
reconnection, which informs us about how fast reconnection can 
take place in a given steady state situation. 

One problem that has been encountered when starting to 
make numerical attempts to investigate the magnetic reconnec- 
tion scenario is that the imposed magnetic resistivity (77) has sig- 
nificant implications for the development of the reconnection re- 
gion. In the case where constant resistivity is imposed, the mag- 
netic structure has a tendency to collapse into a current sheet, but 
it has a difficulty in reaching a steady state situation where there 
is a perfect balance between the inflow and the reconnection pro- 
cess with a constant size current sheet. Instead, the current sheet 
tends to continue growing towards a solution that resembles a 
Sweet-Parker current sheet, with close to "infinite" length and 
close to zero reconnection rate. The way to avoid the infinite 
current sheet is to change the 77 profile. As soon as n is made 
space, time or field dependent, the situation changes and the cur- 
rent sheet reaches a finite extent instead. 

Over the last 10-15 years the focus of investigations has 
changed from 2D reconnection to the more challenging and real- 
istic 3D reconnection processes. It is becoming clear that in 3D 
the possible reconnection scenarios are much more diverse than 
in 2D, and provide a series of new challenges to resolve. The 
process of 3D reconnection is not restricted to regions where the 
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Fig. 1. The structure of a single null. The direction represents 
the spine axis, while the plane defined by e\ and e2 represents 
the fan plane. 



magnetic field vanishes as in 2D, but may also take place in re- 
gions with non-vanishing magnetic field. In this paper we will 
investigate in some detail the reconnection process that can take 
place in the vicinity of a single 3D null point. Null points consti- 
tute critical points of the magnetic field, with the only possible 
class bein g a hyperbolic critical point due to the solenoidal na- 
ture of B fcreenft 19891; lLau & Finniri990t iParnellet al.lll996l) . A 
single pair of field lines approach (recede from) the null from op- 
posite directions, and identify the so-called "spine" of the null. 
A set of field lines radiate away from (approach) the null in a 
surface known as the "fan" of the null (see Fig . [TJ . 

Earlier investigations have shown that reconnection at 3D 
nulls can take place in a number of characteristic manners. First, 
the topological properties of the magnetic flux evolution have 
been shown to be crucially dependent on the orientat ion of the 
electric current at the null dPontin et al.ll2004l 120051) . Second, 
depending on the motions that drive the formation of the cur- 
rent layer at the null, the reconnection may occur either in a 
tube aligned to the spine, a planar layer in the fan, or a layer 
focus ed at the null (Ric kard & Titovl[l996t iPontin & Galsgaard! 
20071). The various reg imes have recently been categorized by 
Priest & Pontinl d2009l) Due to the increased complexity of the 
magnetic structure, analytical solutions to the problem of single 
null reconnection can only be achieved using various simplify- 
ing assumptions, and no-one has at present shown how steady- 
state reconnection takes place at a single null in a compressible 
plasma, or if indeed a steady-state solution is possible. 

In the present investigation a 3D null is stressed by impos- 
ing a systematic driving velocity on the two boundaries through 
which the spine axis of the null penetrates. The driving is of a 
shear type, such that the angle between the spine and fan would 
be expected to ch a nge. U nder the simplistic assumptions made 
by iPriest & Titovl ([1996) this was expected to lead to current 
accumulation in the entire fan plane. Although it appears that 
such a situation would indeed occur in an incompressible plasm a 
(ICraig et al] 1 19951: ICraig & Fabling! Il998l IPontin et al.ll2007bh . 
when plasma compressibility is included, a local collapse of the 
null occurs, destroying the planar nature of the fan surface. A 
localised current sheet forms which is focused around the null 
itself, with the result that the geometry of the magnetic field 



around the null is significantly different from the initial linear 
profile. The resulting mode of reconnection has b een dubbed 
"spine-fan reconnection" by IPriest & Pontinl (120091) (having the 
appearance of a hybrid of early spine and fan models), and in- 
volves magnetic field lines being transported through/around the 
spine and across the fan surface. 

A series of numerical experiments have investigated the 
reconnection process at single null. These have used various 
forms of driving of the system, that ge nerally have perturbed 
the null region for a finite period of time ([Galsgaard et all 2003: 
IPontin et alJ2007blal : IP~ontin & Galsgaardl2007t) . This paper dis- 
cusses the situation where a long time systematic driving is im- 
posed on the system in an attempt to reach a steady state re- 
connection configuration. By continuously advecting the mag- 
netic field in the vicinity of the spine axis, the magnetic field 
is compressed towards the null from two sides. Assuming that 
reconnection takes place, then due to the change in magnetic 
field strength away from the null, one would expect a constantly 
changing amount of magnetic flux to be advected towards the 
null per unit time. As a result, one might expect it to be diffi- 
cult to obtain a true steady-state situation, while a continuously 
evolving quasi-steady-state situation may be more likely. 

In this paper we will discuss, in Section|2] the generic set-up 
of the problem, the numerical approach and the imposed driving 
of the system. In Section [3] the experiments will be discussed 
and the results highlighted. This is followed, in Section [4] by 
a discussion of the consequences of the results for our general 
understanding of the null point reconnection process. Finally, in 
Section[5]we present our conclusions. 

2. Numerical setup 

The magnetic field configuration adopted here is the linear rota- 
tionally symmetric potential magnetic field defined by 

B = B [-2x,y,z\, (1) 

where Bo is a scaling parameter, here of order unity. The null 
is located at (x, y, z) - (0, 0, 0) with the spine axis along the x- 
axis and the fan plane represented by the y — z plane, i.e. x — 0. 
This linear null point is simple to use in numerical experiments, 
but represents a non-physical situation where the magnetic field 
strength continues to increase with distance from the null point. 
In reality non-linear effects will become important at some dis- 
tance and change the linear growth. To investigate the dynami- 
cal evolution of the null, the domain of interest has to be lim- 
ited in size. For this investigation we have limited the domain to 
(x,y,z) = ±(0.5, 1.5, 1.5), and have discretised it using a uniform 
grid with 200x300x300 points. 

The plasma is initially assumed to be at rest and have a con- 
stant density and thermal energy. In the present experiment the 
MHD equations are non-dimensionalised and the density is set 
to unity, while the thermal energy is set to 0.025. This means that 
the surface on which the plasma beta equals unity is an ellipsoid 
with minor axis along the x-axis of length 0.08 and major axis 
perpendicular to this with length of 0.18. Inside this surface the 
plasma beta tends to infinity as the null point is approached. 

The linear nature of the magnetic null point forces a lim- 
itation on the numerical domain size. As a Cartesian domain 
does not represent very well the initial rotational symmetry of 
the magnetic null structure, and because we are going to impose 
a symmetry breaking perturbation, it is a non-trivial problem to 
describe the boundary conditions of the imposed domain as the 
perturbation reaches the boundaries. In the present approach, we 
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Fig. 2. A surface plot of the imposed driving velocity profile. 

force the flow perpendicular to all of the boundaries to be zero, 
and restrict parallel velocities at the boundaries to be non-zero 
only in the two regions where dedicated stressing velocity flows 
are imposed. 

In these experiments we limit the driving of the system to 
two areas on the x -boundaries, by advecting the magnetic field 
in the y-direcfion, with an additional spatial dependence on the 
z-direction. The imposed driver has the following expression, 

V y (y,z) = ±0.5 V d (tanh(Cy - yo)/yh) - tanh((y + y )/y h )) (2) 
x (tanh((z - z )/z h ) - tanh((z + z )/z h )), 

where Vd is the peak driving amplitude, and (y,z) = (±yo, ±Zo) 
represent the locations (at the centre of the hyperbolic tangent 
function) where the driving velocity has half of its peak value. 
The variables y/„ Zh are the half widths of the hyperbolic tan- 
gent functions. This gives a near linear driving velocity between 
(y< z) = (±(yo - yh), ±(zo - Zh)) while the velocity decreases out- 
side this region, asymptotically approaching zero. For these ex- 
periments we use yo = 0.6, yh = 0.2, zo = 0.3 and zi, = 0.2 (see 
Fig.|2]i. This profile is imposed on the two x-boundaries with op- 
posite directions. The driving is switched on at the initiation of 
the experiment and reaches the peak velocity exponentially with 
a given time scale. 

The experiments are conducted usi ng the 3D non-ideal MHD 
code bv lNordlund & Galsgaardl(ll997l) . This is a high order finite 
difference code using staggered grids to maintain conservation 
of physical quantities. The interpolation operators are fifth order 
in space while the derivative operators are sixth order. The so- 
lution is advanced in time using a third order explicit predictor- 
corrector method. Here we intend to derive scaling relations, so 
we adopt a constant rj model. 

3. Experiments 

The aim of the experiments is to investigate the possibility of ob- 
taining steady state reconnection at a 3D null point. We also seek 
to derive scaling relations to enable the extension of the results 
to higher Reynolds numbers. In what follows we investigate the 



Table 1. Experiments investigating the impact of the imposed 
boundary driving velocity, Vd, and the plasma resistivity, tj. 



Name 




n 




Peak J b 


Peak V c 


(1 

II 


Al 


1 




0.100 


8.0 


0.75 




A2 


1 


icr 3 


0.100 


41.6 


0.83 


0.033 


A3 


1 


1(T 4 


0.100 


140.0 


0.93 




B 


1 


10" 3 


0.075 


38.4 


0.73 


0.026 


CI 


1 


lO 2 


0.050 


5.0 


0.34 


0.040 


C2 


5 


io- 3 


0.050 


9.2 


0.42 


0.034 


C3 


1 


1(T 3 


0.050 


33.0 


0.58 


0.023 


C4 


5 


10 4 


0.050 


52.0 


0.68 


0.019 


C5 


1 


ltr 4 


0.050 




0.89 


0.0058 


C6 







0.050 








D 


1 


ltr 3 


0.020 


19.3 


0.28 


0.013 


E 


1 


ltr 3 


0.005 


3.5 


0.082 


0.0037 



a imposed boundary driving velocity 

b peak current value 

c peak outflow velocity 

d peak integrated parallel electric field 
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Fig. 3. Shearing of a 2D X-point. Left: sheet length as a function 
of time for domain [x, y] e [+0.5, +1.5] (dashed line) and [x,y] e 
[±1, ±3] (solid line). Right: peak current as a function of time for 
the same two domain sizes. 

effect of varying two different global parameters of the system, 
namely the magnetic resistivity (77) and the imposed driving ve- 
locity (V d ). 

An overview of the various simulation runs is given in 
Table Q] together with three characteristic parameters for the ex- 
periments that will be discussed later. Adopting characteristic 
parameters for the experiment, we can determine the magnetic 
Reynolds numbers for these experiments Rm = V a Ljr\, where 
V a = 0.5 is the Alfven velocity at the spine axis at the top/bottom 
boundary, L — 1 is the distance between the x-boundaries and rj 
is given in Table Q] This leads to values for Rm between 50 and 
5000. 

The following subsections discuss different aspects of the ex- 
periments. We begin with a general description of the dynamical 
evolution, first discussing a 2D setup as a reference for the 3D 
experiment, since our knowledge of the behaviour in the 2D case 
is well established (at least within the MHD framework). 

3.1. Evolution of a 2D system 

We first briefly describe the evolution of a sheared 2D null 
point. The magnetic field is taken to be B = [—x,y, 0] and 
the z-dependence of the driver is suppressed. As the driving is 
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switched on two MHD waves propagate towards the null point 
from either of the driving boundaries. Close to the null they col- 
lide and a strong current layer forms in the region where the two 
separatrices collapse very close to one another. The length of the 
current sheet and the maximum current density (and therefore 
in 2D the reconnection rate) gradually build to peak values, be- 
fore beginning a slow decay (see Fig . [3j ■ During this period, the 
configuration can be described as being in a quasi-steady state. 

There are a number of key points to note. First, the halt 
in the growth of the current layer (both length and modulus), 
is due to the finite size of our computational domain. If we 
repeat the simulations in a larger computational domain, with 
an extended driving region, the current layer continues to grow 
for substantially longer, before again being constrained by the 
boundaries (compare the full and dashed lines in Fig. |3j. This 
is a well-known effect in 2D - that for uniform resistivity a 
'Sweet-Parker' type current layer develops, which continually 
grows as the forcing continues. The limitation on the current 
layer in these simulations is therefore equivalent to the limitation 
imposed when reconnection i s forced in an in itial Harris sheet 
with periodic boundaries (e.g.lBirn et al. 2005). Another impor- 
tant point to note is that, were the driving really constant in our 
simulations, we would not expect such a noticeable decay of the 
current layer length and modulus after the peak is reached. The 
decay is due to a reduction in the quantity of flux being forced 
into the sheet per unit time as the simulations progress. This in 
turn is due to the evacuation of flux from the region upstream of 
the inflow, as flux from that region is constantly fed through the 
sheet. The above properties all have analogues in the 3D system, 
as discussed below. 
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Fig. 4. Current amplitude at a number of characteristic instances 
in the dynamical evolution. The individual frames are scaled to 
their local dynamical range to enhance the visual appearance. 
The images represent the z — plane for run "C3" at times t 
(0,4. 1.0.2.2,4,4.9.2. 19.0). 



3.2. Qualitative evolution of the 3D system 

Initiation of the driving on the two x-boundaries launches two 
wavefronts that propagate through the magnetised plasma. As 
the driving speed is slowly ramped up the wavefronts ini- 
tially propagate towards the fan plane as simple plane Alfven 
waves (with their x speed being independent of y and z). With 
the driving continuously stressing the magnetic field, the field 
lines behind the wavefront are systematically advected in the y- 
direction. This breaks the initial symmetry of the field, gener- 
ating an asymmetry of the field on either side of the advected 
spine line. This causes the wavefront to tilt as its speed becomes 
dependent on y, due to an increase in wave speed in the region 
where the field lines are stretched, Fig. |4] 

From the initial phase of the evolution, two points should be 
noted. First, the width of the current front decreases while its 
intensity increases on its way towards the fan plane. This is an 
effect of the decreasing Alfven velocity in the x-direction, that 
for a pure Alfven wave implies that it can never reach the fan 
plane of the null point in a strict mathematical sense. This results 
in a pileup producing a strong current front. As this occurs, the 
wavefront also spreads out increasingly fast along the y- and z- 
directions as it approaches the fan plane, due to the field geome- 
try. One would expect this to act to reduce the current intensity in 
the front, were this effect the dominant one. However, due to the 
compression of the plasma, the wave is not a pure Alfven wave, 
but couples to both slow and fast mode waves. The perturbation 
can therefore propagate across the magnetic field lines and ad- 
vance towards the null point. This type of behaviour was seen in 
Galsgaar d et al.l (120031) . where the null point was perturbed by 
two rotational twist waves centred on the spine axis. Here the 
fast mode wave wraps around the null and concentrates part of 



the en ergy of the perturbation there (see also iMcLaughlin et al.l 
120091 who studied the 2D case). 

As the evolution proceeds, the two wavefronts initiated at op- 
posite x-boundaries approach each other and eventually collide, 
with the two non-rotationally-symmetric current concentrations 
compressing into one central current distribution centred on the 
null. Within the resulting current layer, the spine and fan have 
collapsed towards one another, and the current spreads along the 
fan both along and across the driving direction. The tilt and mag- 
nitude of the arising current concentration depend on the im- 
posed driving speed and imposed 77 as will be discussed later. 
Typically the tilt of the current sheet overshoots initially and the 
sheet experiences a few oscillations before settling to a stable 
angle. Due to the imposed constant 77, reconnection takes place 
in the current sheet, and a near steady-state balance between the 
flux advected into the current sheet and the reconnected flux ex- 
iting the sheet through a reconnection jet is reached (see below). 
This picture is maintained with only minor changes until the ter- 
mination of the experiments. 

3.3. Steady state evolution 

From the 2D experiment it is seen that the evolution reaches a 
quasi-steady-state. In order for it to be physically meaningful to 
search for scaling relations, we must first check if such a quasi- 
steady-state also exists in the 3D simulations. To determine this 
we examine some of the charac t eristic para meters addressed in 
the Sweet-Parker (lParkerlll957t ISweetll 19581) 2D reconnection 
model. 

In order to define the characteristic size of the current layer 
we define the boundary of the layer to be the surface at which 
the current modulus |J| falls to 50% of its peak value (this peak 
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Fig. 5. The time evolution of the magnetic field in the inflow re- 
gion (top left), the plasma inflow velocity (top right), the outflow 
velocity relative to the Alfven velocity in the inflow region and 
the ratio between the inflow (lower left) and outflow velocities 
and sheet dimensions defined by Lvj„/(lv ol „) (lower right). For 
the run C3. 
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time time 



Fig. 6. The four panels show the time evolution of the peak value 
of the current in the sheet (top left), the peak velocity of the 
reconnection jet (top right), the length of the current sheet (lower 
left) and finally its angle relative to the y-axis (lower right), for 
the run C3. 



occurring at the null). In 3D this provides three characteristic 
length scales that are here listed in descending order; The length 
(L) in the plane of the shear perturbation (here the xy-plane), the 
extent of the sheet in the direction perpendicular to this plane 
and parallel to the current at the null (here the z direction) which 



we denote here by the width (w), and finally the thickness, I, of 
the sheet (measured in the z = plane). 

First we estimate the amount of magnetic flux being trans- 
ported into the sheet - for which one needs to know both the in- 
flow velocity perpendicular to the current sheet and the magnetic 
vector aligned with the sheet. For the purposes of making such 
an estimate, we measure these values at the centre of the inflow 
region just outside the current sheet. We also consider the out- 
flow jet velocity, measured toward the outflow edge of the cur- 
rent sheet. Fig. [5]shows four combinations of these variables as 
functions of time from experiment "C3". The top left frame rep- 
resents the component of the magnetic field parallel to the cur- 
rent sheet. It shows an initial build-up followed by a near linear 
decline in magnitude. The latter is a consequence of the deple- 
tion of the magnetic field in the inflow region due to the limited 
length scale of the driving in the y-direction. The top right panel 
shows the inflow velocity, which is found to reach a stable level 
with only a slow decay with time. The lower left frame shows 
the ratio of the outflow velocity and the Alfven speed in the in- 
flow region (akin to the reconnection rate in the Sweet-Parker 
model), which reaches a nearly stable value at the same time as 
the magnetic field component reached its peak value. Finally the 
lower right frame shows the ratio of the length scale of the inflow 
region times the inflow velocity and the thickness of the outflow 
region times the outflow velocity. In the 2D Sweet Parker anal- 
ysis this relation constitutes mass conservation, and here we can 
see that it reaches a stable value declining only slowly with time. 
It is interesting to note that this stable value is somewhat greater 
than 1 (the 2D value), being closer to a value around 1.5. This is 
likely due to the fact that the outflowing plasma may access the 
third dimension in 3D, rather than being constrained to flow out 
in the plane of the inflow. 

In Fig.|6]the top left panel shows the time development of the 
maximum current density in the domain, which reaches a peak 
around t = 7 followed by a slow decay. The top right panel repre- 
sents the time development of the jet velocity, which is slightly 
different from the current, taking longer to reach its maximum 
value. This is expected since the velocity is a collective effect 
of the Lorentz tension force and gas pressure in the newly re- 
connected field lines, acting on a "dense" plasma and therefore 
needing some additional time to accelerate the plasma to it cruis- 
ing velocity. The lower left panel shows how the sheet length is 
very stable in time, while the lower right panel represents the tilt 
angle of the sheet that shows a slow decreases with time. 

From the above discussion it appears that the effect of the 
continual driving is to force the null structure into a near steady 
state, in which the reconnection region maintains its character- 
istics for a period of time much longer than the Alfven travel 
time along the current sheet (Fig. [5] and Fig. [6}. The qualitative 
behaviour discussed above is also observed in the other experi- 
ments listed in Table [T] above. We therefore go on in the follow- 
ing section to derive scaling relations between various parame- 
ters in the simulations. 

3.4. Scaling relations 
3.4.1 . Peak current 

It is well known that the peak current J ma x reached in a numer- 
ical experiment depends on two key parameters, one being the 
numerical resolution. If the imposed 77 is too low then the nu- 
merical diffusion dominates and will limit J nuu . The other key 
parameter is the chosen value of the resistivity 77, which limits 
the current growth through the "correct" evolution of the induc- 
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Fig. 7. Scaling relation between the imposed 77 and the thickness 
of the current sheet (measured in units of the grid spacing in the 
jt-direction), for the runs CI - C6. The full line represents sheet 
thickness oc 77 035 . 
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Fig. 8. The scaling relation of the peak current with the im- 
posed 77 value. The data points represent two series of experi- 
ments with different imposed driving velocities. The diamonds 
and the associated full line represent data with V £ /=0.05. The 
triangle represents the corrected 77 for the left most point as in- 
dicated in the previous graph, with the dot-dashed line being the 
correspondingly altered scaling. The crosses and the dashed line 
represent the data with V c /-0.l (with the star showing the cor- 
rected 77 value). The dashed line represents Eq. 3, the full line 
Eq. 4 and the dot-dashed line Eq. 5. 



tion equation. With sufficient numerical resolution one can re- 
solve the structure of the current sheet and it becomes possible 
to make a meaningful analysis of the impact of changing 77. In 
this investigation we have run simulations with a limited selec- 
tion of 77 values. 

One way to investigate whether the imposed resistivity dom- 
inates over the numerical resistivity is to measure the width of 
the current sheet as a function of the imposed 77 value. In the 
incompressible case iHeerikhuisen & Craig (2004) showed that 
this gives rise to a power law dependence with a power of 0.5 
(although we note that the morphology of the current sheet is 
somewhat different in their incompressible models to what is ob- 
served here). Fig. [7] shows the relation for the "C" experiments 
given in Table Q] The full line is a power law, but with a power 
of 0.35 for the higher 77 values. In the plot the diamonds repre- 



Fig. 9. The scaling relation between the imposed driving veloc- 
ity and the peak current for the experiments with 77 = 10~ 3 listed 
in Tabled 

sent the imposed values of 77, with the point furthest to the left 
representing a case where 77 was set to zero. This indicates that 
the code introduces an efficient 77 when the width of the sheet be- 
comes on the order of the numerical stencil. If one assumes that 
the power law dependence can be extended, then an asymptotic 
value of the numerical 77 can be estimated for the given reso- 
lution. Correcting the values of 77 correspondingly for the runs 
C5 and C6 (two left most points on the graph) brings those data 
points onto the line representing the scaling relation (diamonds 
in the plot). In the following we will investigate if this correc- 
tion of the imposed resistivity will have any implication on the 
scaling relations. 

When data from the experiments listed in Table Q] are used, 
scaling relations of the peak current with 77 can be obtained. 
Fig.[8]shows plots of peak current versus 77 for two different driv- 
ing velocities. The triangle represents the "corrected" 77 value for 
the run "C4" and the star is the corrected value for "A3". The 
scalings, corresponding to the lines plotted in the Figure, are 

Jmca ■. = 0.46 77- 062 for 7^ = 0.1, (3) 
J max = 0.14 77~° 78 for V d = 0.05 - without 77 correction{4) 
J max = 0.1 77-° 85 for V d = 0.05 - with 77 correction. (5) 

Here one can see two different comparable fits for the slow 
(Vd = 0.05) driving experiments, with the scaling based on the 
corrected value leading to a slightly stronger dependence of J max 
on 77. However, the correction of the 77 value for the faster driver 
does not seem to lead to a linear scaling relation. 

We now turn to the scaling of the peak current with the driv- 
ing velocity, for a constant 77 value, which is investigated us- 
ing the 77 = 10~ 3 experiments listed in Table Q] The result is 
shown in Fig. [9] where it is seen that while the peak current 
clearly depends on the imposed driving velocity, there is no sim- 
ple (e.g. linear or exponential) relation that can be obtained to fit 
the data. 

3.4.2. Current layer dimensions and jet velocity 

The dimensions and tilt angle of the current sheet have been 
measured for the same experiments. They are found to be inde- 
pendent of 77 for a constant driving velocity, to within our mea- 
surement accuracy. On the other hand, varying the driving veloc- 
ity while maintaining a constant value for 77, we find that there 
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Fig. 10. The top frame shows the time dependent evolution of 
the peak velocity, for the runs CI (dot-dashed), C2 (dashed), 
C3 (triple-dot-dashed), C4 (long-dashed) and C5 (solid) (see 
Table[T]i. The middle frame shows the peak velocity as a function 
of the imposed driving velocity. The lower frame shows the peak 
velocity as a function of the imposed resistivity, 77. The full line 
represents the scaling using the original values, Eq.(6), while the 
dot-dashed line included the corrected 77 values, Eq.(7). 



are significant changes of these dimensions, indicating that they 
are largely dependent on the amount of imposed stress. However, 
as discussed above, the present domain size and the extent of the 
driving region limit the current sheet from evolving freely, and 
thus with the present setup we are unable to provide scaling re- 
lations for these quantities. 



The top frame of Fig.[10]shows the maximum outflow veloc- 
ity as a function of time for different driving velocities Vd- This 
shows a significant variation both with time and imposed driving 
speed. The middle frame shows the peak outflow jet velocity (for 
all time) as a function of the imposed driving velocity. Here we 
can see that we have a non-linear dependence, which is also not 
fit well with an exponential scaling, but follows a similar pattern 
to the scaling of J max with Vd- Finally, the lower frame of Fig.fTOl 
shows the jet velocity as a function of the imposed resistivity, 
where we find approximately power law scalings, with 

Vpeak - FT 7 ?" 22 _ without 77 correction, (6) 
8.1 

V pea k = g~3 7 7"°' 25 ~~ with 1 correction. (7) 

The first scaling is for the case where we use the manually ap- 
plied 77 values for the experiments, while in the second the cor- 
rected values are included. 

3.5. Reconnection rate 

From established theory (ISchindler et al.l 119881 iPriest et al] 
2003) we know that the correct measurement of the reconnec- 
tion rate in 3D is given by the maximal value of <t>, the inte- 
grated electric field parallel to a magnetic field line (E\\), where 
this maximum is taken over all field lines which thread the non- 
ideal region (assumed to be spatially localised). To determine 
the reconnection rate we therefore perform a systematic field line 
tracing for a collection of field lines that thread the current sheet, 
and seek the maximum in the integrated parallel electric field. In 
practice this occurs along a field line that is located very close 
to the z-axis (it would be expected to lie exactly alon g the z-axis 
based on the exactly symmetric analytical model of lPontin et al.l 
(2005)). In these experiments a constant 77 has been used, which 
implies that E\\ — tjJ\\, so that 

<D = J (77J • B)/|B|<77, (8) 

where we integrate along the given field line starting close (e) 
to the null out to a finite distance, L. To get the correct recon- 
nection rate the integral has to be made in th e z-direction in both 
senses, to ±L. In previous simulations (e.g. iPontin et al.l l2007a) 
this was a relatively simple process, as the diffusion region was 
limited to a finite size. However, the situation here is different. 
Due to the continual driving, the current not only has a strong 
component close to the null, but there also exists an extended 
(though weaker) concentration of current in the fan plane that 
reaches all the way to the z-boundaries. This is demonstrated in 
the top image in Fig. QT| which shows a ID plot of the paral- 
lel electric field along the z-axis (which is not a magnetic field 
line). In the lower frame a 2D (y, z) representation of the same 
quantity is presented, which again shows the extended parallel 
electric field. Plotting the same quantity in a plane perpendicular 
to this shows that the parallel electric field is by contrast limited 
to a very narrow layer in the jc-direction (compare with Fig. 3). 
The actual reconnection process is therefore not limited to the 
local region around the strong current sheet, but extends all the 
way to the imposed boundary. This clearly implies that the value 
obtained using Eq. © depends on the choice of L. For the results 
below the value of L has been limited to L = 1 .4, such that the 
domain boundary is not reached. 

Performing the integration defined in Eq. ((8) on the various 
data sets, we obtain the results shown in Fig. [12] The left frame 
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Fig. 11. The top panel shows the parallel electric field (E\\ = 
E ■ B/|B|) along the z-axis for the run "C3". The bottom panel 
shows a 2D image of the parallel electric field in the x — plane. 

shows the peak value of as a function of rj. If we correct the re- 
sult obtained with the lowest rj value (run C5, left most diamond 
in the plot) relative to the results obtained from Fig. [7] this point 
is shifted to the location marked by the triangle in the plot, that 
falls nicely on the straight line scaling relation that fits the other 
data points. This simple correction seems to work here, most 
likely because the current sheet in the diffusion region is all the 
way just being resolved with a minimum number of gridpoints. 

In the right hand frame, an indication of a scaling relation 
between the driving velocity and reconnection rate can be seen, 
although with some fluctuations around the obtained line. The 
two power law relations are given by: 

<D = J~£||<« = 0.126 77 025 for V d = 0.05 (9) 
<D = J E {l dl = 0.15 V% 65 for7/=10- 3 (10) 

4. Discussion 

The experiments described above were conducted to investigate 
whether a steady-state magnetic reconnection process could be 




0.0001 0.0010 0.0100 0.01 0.10 

eta v drive 



Fig. 12. The left panel shows the dependence of the recon- 
nection rate on the imposed rj with a constant driving velocity 
(Yd = 0.05), with the triangle showing the corrected value for 
the run "C5" as discussed above, Eq.(9). The right panel shows 
the dependence on the reconnection rate on the imposed driving 
velocity, for constant 77 (77 = 10~ 3 ), Eq.(10). 



set up in 3D by continually stressing a symmetric null point. The 
results presented in Section l3~3l show that a quasi-steady state ex- 
ists where the inflow to the diffusion region is balanced by the 
ongoing magnetic reconnection (this is true for all values of the 
imposed driving velocity and rf). Due to the magnetic geometry, 
the inflow magnetic field and the plasma inflow and outflow ve- 
locities (v,„ and v out ) are not constant during this quasi-steady 
evolution. However, ratios between characteristic quantities fa- 
miliar from 2D models do settle to approximately constant val- 
ues. For example, the ratio of v out to the inflow Alfven speed 
is approximately fixed. So is the ratio Lvi„/lv out , which charac- 
terises mass conservation in the sheet and settles to a constant 
value of around 1.5 (rather than 1 as in 2D, since the outflowing 
plasma may access the third dimension rather than being con- 
strained to flow out in the plane of the inflow). It should be noted 
that in each of the simulations the quasi-steady configuration 
obtained is strongly influenced by the restricted dimensions of 
both the numerical domain and the region on which the bound- 
ary driving is imposed. If it was possible to extend the domain 
and driving region indefinitely, indications are that the current 
layer would extend indefinitely, until reaching a threshold for 
some instability such as the tearing instability. 

We have investigated the dependence of different character- 
istic quantities in the quasi-steady state on both the driving ve- 
locity Vd and imposed resistivity 77. Considering first the depen- 
dence on Vd, we found that the current maximum and outflow 
velocity increase as Vd increases, although not with a simple 
scaling relation. It is natural that these quantities increase when 
the rate at which flux is forced into the sheet is increased, and 
also that they behave in the same way since the outflow jet is 
accelerated by the J x B force. 

For the simulation runs with different values of 77 we have 
been able to determine approximate scaling relationships. The 
peak current appears to scale like J max ~ rj" 6 - 7]~ & , which 
is consistent with the rf^^ s caling expected for pla nar incom- 
pressible fan current layers dCraig & Fabling|[l998l) . The peak 
outflow velocity on the other hand follows an approximate scal- 
ing Vp e ak ~ i]~ 25 , which mirrors the scaling of the reconnec- 
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tion rate, as we discuss below. It is also clear from Fig. |7]that 
the width of the current sheet depends on the imposed 77 value, 
and the graph indicates that as 77 approaches zero the sheet 
reaches a minimum thickness, controlled by numerical diffu- 
sion. Results for large 77 indicate a power law dependence for 
the thickness as a function of 77, though with a smaller expo- 
nent than found in the incompress ible models (Craig & Fabling 
1998; Heerikhuisen & Craig 2004). One can use this to estimate 
the value of the numerical diffusivity, which takes effect only 
when structures collapse to the resolution limit, while its value 
decreases rapidly as the structures increase in size. Making use 
of this information to correct the scaling parameters in different 
plots is not simple, as indicated by the additional points included 
in a number of figures. In some cases it seems to provide a scal- 
ing relation, while in other cases, it makes no coherent correc- 
tion to the result. The only way to provide better information for 
these values of 77 is through experiments with even higher nu- 
merical resolution. However, doubling the numerical resolution 
locally makes the experiments much more expensive. To make a 
serious difference in the scaling at least a factor of 10 variation 
in resolution is required - which gives a factor of 10 4 in required 
computing time. It is therefore a significant computational chal- 
lenge to improve these scaling relations significantly for lower 77 
values. 

From 2D reconnection theory the simple Sweet-Parker 
model gives a scaling of the reconnection rate that depends on 
77 5 . With realistic parameters for astrophysics plasmas this re- 
connection rate is far too slow to make it an attractive mechanism 
for magnetic energy conversion, in for instance solar flares. It 
has therefore since then been a task to find reconnection setups 
that depends less strongly on 7/ and therefore evolve much faster. 
Taking the scaling rate for the reconnection rate given in Eq. (0, 
the driven fan-spine reconnection is found to scale with a smaller 
power of 77, namely 0.25. This scaling for the reconnection rate 
is significantly faster than the slow Sweet-Parker scaling, though 
slower than the fast Petcheck reconnection rate that scales with 
the logarithm of tj. However, it is worth pointing out that over our 
range of 77 covering two orders of magnitude (77 = 10~ 2 - 10~ 4 ) 
it is difficult to distinguish the power law and ln(?7) scalings. 

The obtained scaling relations should be compared with 
those obtained for the same setup, but with an impulsive 
(i.e. temp orally localis e d) pe rturbatio n of the null point, de- 
scribed by iPontin et al.1 (l2007al) and iPriest & Pontinl (120091) . In 
these investigations, the boundary driving was applied for a fixed 
period of time and then reduced to zero. It is worth noting that 
this implies that the net displacement of the spine - under the 
assumption of an ideal evolution - increases as the driving speed 
increases. Their analysis showed that for fixed 77, the peak current 
and reconnection rate increase linearly with the driving velocity, 
while the length and width of the sheet (L and w above) actually 
decrease linearly. This seems counter-intuitive, but is a result of 
the fact that the collapse of the field in the vicinity of the null is 
more efficient for stronger driving. Thus, the scaling of the cur- 
rent and reconnection rate follow similar patterns to those seen 
above in the continually driven case, however the sheet dimen- 
sions scale in the opposite way. Turning now to the scaling with 
resistivity 7/ (for fixed driving speed) in the impulsively driven 
case, it was found that as 77 decreases, the peak current increases 
and the reconnection rate decreases, with a rate between power 
law and logarithmic, while here we may argue in favor for a 
power law dependence. Clearly there are some similarities, but 
also some differences, between the two different driving config- 
urations. 



5. Conclusions 

The experiments have shown that it is possible for a 3D null 
point to reach a quasi-steady state, where the characteristic pa- 
rameters of the current sheet are controlled by a combination of 
the imposed driving velocity and the effective constant resistiv- 
ity. It is likely that other parameters, such as the dimensions of 
the plasma-/3 = 1 surface also play a role. It is found that scaling 
relations between characteristic parameters do not always follow 
simple exponential or power law expressions over the full range 
of rj values. This is partly due to the limited numerical resolution 
in the simulations, which means that the code imposes its own 
diffusion when structures becomes too small. Only with substan- 
tially higher numerical resolution will it be possible to push the 
reliable range of 77 to significantly smaller values. However, this 
preliminary study indicates that the reconnection rate scales rel- 
atively weakly with the plasma resistivity, as 77 1 / 4 . 

The scaling relations found here are different from the 
ones derived for the impulsive perturbation of the single null. 
However, we find that even changing the amplitude of the fixed 
driving velocity changes the scaling with resistivity. It is there- 
fore not easy to make any predictions regarding the global be- 
haviour of 3D nulls in realistic situations: too many parameters 
seem to influence the dynamical evolution, with the scalings be- 
ing linked in ways that appear to be non-linear, and that we are 
therefore yet to unravel. Further to this, test experiments indicate 
that some of the obtained parameters are likely to change when 
the size of the imposed driver is changed. 

It is clear from the results that more effort has to be invested 
in understanding null reconnection and its scaling behaviour be- 
fore we are able to make any quantitative predictions regarding 
processes occurring at 3D nulls in the Solar atmosphere or the 
Earth's magnetosphere. 
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